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FINAL  REPORT  OF  NRL  RESEARCH  IN  SUPPORT  OF 
DOE  CHARGED-PARTICLE  INERTIAL-CONFINEMENT-FUSION 
PROGRAM  FOR  FY  77-78 


I.  INTRODUCTION  TO  THE  RESEARCH  PROGRAM 

During  the  reporting  period.  Plasma  Physics  Division  research  has  contributed  to  the 
Department  of  Energy  charged-particle-beam  inertial  confinement  fusion  program  in  two  major 
areas: 


•  electron-beam  diode  physics  investigating  the  achievement  of  predictable  and 
reproducible  intensely-pinched  beams 

•  charged-particle-beam  transport  physics  for  pellet-ignition  systems  employing  elec¬ 
tron,  light-ion  and  heavy-ion  beam  drivers. 

This  report  documents  research  in  these  areas  which  has  not  appeared  in  print  in  Sec.  1I-V. 
Published  papers  and  reports  which  support  the  FY  77-78  DOE  program  are  listed  in  Sec.  VI 
and  are  referenced  as  subsection  letter  headings. 

A.  ELECTRON-BEAM  DIODE  RESEARCH 

Electron-beam  diode  physics  research  carried  out  during  the  reporting  period  can  be 
divided  into  two  categories: 

•  Detailed  experimental  and  theoretical  investigations  of  the  behavior  of  conven¬ 
tional  tapered,  pinched-electron-beam  diodes. 


Manuscript  submitted  April  29.  1980. 


MOSHER  et  al. 


•  Preliminary  experimental  and  theoretical  investigations  of  a  novel  diode 
configuration. 

Research  in  the  first  category  provides  scaling  relations  for  the  electrical  characteristics  of 
diodes  and  documents  electron-electrode  interaction  information  required  to  achieve  highly- 
focused  beams.  Work  in  the  second  category  documents  the  behavior  of  a  diode  configuration 
with  improved  coupling  characteristics  to  low-impedance  generators  and  high  current  density 
capability. 

Papers  VI  A-C  cover  the  first  category.  The  time-dependent  impedance  behavior  of 
large-aspect-ratio  pinched  electron-beam  diodes  with  no  external  magnetic  field  is  investigated 
using  the  Gamble  I  and  Gamble  II  generators  in  Sec.  VI  A.  The  experimental  results  are  com¬ 
pared  to  predictions  of  a  time-dependent  computer  code  which  incorporates  several  important 
mechanisms  that  affect  diode  impedance.  A  discussion  of  these  mechanisms  including  experi¬ 
mental  evidence  of  their  existence  and  their  effects  on  the  observed  diode  impedance  is 
presented.  A  simple  model  for  the  current  is  used  in  a  time-dependent  analysis  of  the  data  in 
order  to  demonstrate  qualitative  effects.  The  rise  time  of  the  electric  field,  the  cathode-plasma 
velocity,  the  time-dependent  generation  of  ion  sources  from  the  anode,  and  the  special  diode 
geometries  (tapered  cathodes)  are  shown  to  affect  the  diode  impedance  in  good  agreement 
with  theory. 

The  effect  of  different  anode  materials  on  the  size  of  very  intense  pinches  is  studied  using 
the  radial  structure  of  the  bremsstrahlung  from  a  tightly  pinched  relativistic  electron  beam  in 
Sec.  VI  B.  Time-dependent  data  are  obtained  which  point  out  the  influence  of  scattering  by  the 
anode  which  electrons  experience  near  the  pinch  axis.  The  scattering  broadens  the 
bremsstrahlung  FWHM  for  pinches  with  high-Z  (80)  anodes  to  6  mm,  compared  to  4  mm  for 
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low-Z  (13)  anodes.  When  anodes  composed  of  low-  and  high-Z  elements  in  different  annular 
regions  are  used,  translation  of  backscattered  electrons  enhances  the  current  density  upon  the 
low-Z  regions  and  decreases  it  upon  the  high-Z  regions. 

In  paper  VI  C,  the  EUV  emission  from  10-  to  100-mg/cm2-thick  planar  targets  (Z  —  6- 
79)  exposed  to  60-nsec  duration  focused  relativistic  electron  beams  was  examined  at  current 
densities  of  0.4  -  1.0  MA/cm2.  Spatially  resolved  spectra  were  obtained  with  a  stigmatic  near- 
normal  incidence  spectrograph,  in  the  wavelength  range  500-2100  A,  along  axial  radial  direc¬ 
tions.  Time-resolved  EUV  data  were  obtained  with  windowless  photoelectric  diodes  operated 
unliltered  and  filtered  with  0.5-^m-thick  aluminum.  The  axial  views  show  copious  emission 
from  the  diode  gap  associated  with  late-time  impedance  collapse.  Radial  views  show  continuum 
radiation  which  is  emitted  from  the  region  of  the  target  exposed  to  the  focused  REB  and  line 
radiation  which  originates  from  a  lower-density  (approximately  10”  cm-3)  more-extended  late¬ 
time  plasma  resulting  from  foil  expansion.  The  region  of  the  foil  exposed  to  the  focused  REB 
radiates  as  a  blackbody  source  with  a  peak  temperature  of  2-6  eV.  These  temperatures  are  a 
factor  of  2  below  the  calculated  temperatures  for  classical  deposition  in  the  foils  exposed  to  the 
focused  beam  assuming  no  hydrodynamic  foil  expansion. 

A  radial-injection  electron  beam  diode  configuration  is  discussed  in  Sec.  II.  It  combines  a 
bias-current  pinch  and  an  ion-induced  pinch.  This  technique  allows  one  either  to  minimize  the 
ion  current  or  to  concentrate  it  near  the  axis.  The  radial  injection  was  attempted  as  a  smooth 
transition  for  electrons  from  the  cathode  to  equipotential  surfaces.  A  comparison  of  analytic 
theory,  numerical  simulation  results,  and  experiments  indicates  that  bias-current  pinching  via 
paravector-potential  (PVP)  flow  was  realized  without  reduction  of  ion  flow  in  the  large  area 
regions.  Detailed  geometries  of  diodes  specifically  designed  to  allow  a  smooth  transition  from  a 
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PVP  flow  at  large  radii  to  an  ion  induced  pinch  flow  at  small  radii  were  studied.  Electron  pinch¬ 
ing  from  1 10  mm  to  3  mm  was  observed  with  corresponding  three  order  of  magnitude  current 
density  enhancement. 

B.  CHARGED-PARTICLE-BEAM  TRANSPORT  RESEARCH 

During  the  reporting  period,  research  for  electron,  light  ion  and  heavy-ion  beams  was  car¬ 
ried  out  in  order  to  determine  the  feasibility  of  transporting  focused  or  focusing  beams  through 
a  pellet-fusion  reactor  chamber.  The  investigations  carried  out  include: 

•  High-current  electron  beam  transport  experiments  on  the  Gamble  I  and  Gamble 
II  generators. 

•  Theoretical  investigations  of  focusing  and  transport  of  light-ion-beams  in  external 
magnetic  fields. 

•  Experiments,  theoretical  and  computational  analyses  of  transport  of  light-ion- 
beams  in  channels,  the  MHD  response  of  the  channel,  and  microinstabilities  of 
the  beam-channel  system. 

•  Theoretical  assessment  of  heavy-ion-beam  final  transport  with  respect  to  beam 
neutralization  and  microinstabilities. 

•  Prediction  of  and  estimates  for  enhancement  in  ion  stopping  power  in  ionized 
media. 

In  Section  III,  the  propagation  physics  of  hot  electron  beams  in  plasma  discharges  driven 
by  exploding  wires  is  described.  This  extended  previous  work  carried  out  at  Sandia  Labora¬ 
tories  into  a  regime  of  higher  beam  currents  (up  to  300  kA  injected)  and  investigated  the  effect 
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of  channel  current  and  diode  voltage  on  beam  propagation.  The  effects  of  scattering  and  axial 
electric  fields  on  electron  orbits  as  well  as  the  limitations  on  the  maximum  beam  current  and 
beam  current  density  that  can  be  propagated  was  investigated.  The  experiment  demonstrated 
that  high-power-density  electron  beams  may  be  efficiently  transported  using  current-carrying 
air-plasma  channels.  It  was  shown  that  certain  restrictions  apply  to  their  use,  the  most  impor¬ 
tant  being  limitation  of  the  current  density.  This  restriction  may  be  overcome  to  a  certain 
extent  by  increasing  the  channel  temperature.  This  might  be  accomplished  by  using  a  lighter 
background  gas. 

Research  on  transport  and  focusing  of  light-ion  beams  in  external  magnetic  fields  is  sum¬ 
marized  in  Sections  IV,  VI  D,  and  VI  E.  In  the  reflexive  focusing  scheme  of  Secs.  IV  and  VI 
E,  focusing  on  the  pellet  is  achieved  by  employing  a  weak  solenoidal  field  in  the  transport  sec¬ 
tion,  and  introducing  the  beam  at  small  angles  to  the  field.  The  large  gyro-radius  motion  of  the 
ions  then  brings  the  beam  to  focus  meters  downstream.  A  drawback  of  this  technique  is  that 
the  location  of  the  focus  moves  if  the  ion  energy  varies  in  time.  Thus,  time-of-flight  bunching 
of  the  beam  cannot  be  achieved. 

In  Section  VI  D,  magnetic  compression  of  a  cold,  charge-neutralized  ion  beam  due  to  pas¬ 
sage  through  a  slowly  varying  magnetic-mirror  field  is  studied.  Conditions  on  the  field  and 
beam  parameters  are  derived  which  permit  the  beam  to  be  compressed  to  order-of-magnitude 
higher  current  densities  while  preserving  the  laminar  nature  of  the  flow.  Effects  associated  with 
rotation,  beam  convergence,  and  velocity  shear  in  the  incident  beam,  and  those  due  to  self¬ 
fields  are  included  in  the  analysis.  It  is  shown  that  this  approach  is  of  limited  value  for  focusing 
to  pellet  dimensions  because  of  self-field  effects  and  the  need  for  an  unrealisticly  large  applied 
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The  channel  transport  technique  used  for  electron  beams  is  extended  to  light  ion  beams  in 
Section  VI  E.  By  confining  a  beam  already  focused  to  pellet  dimensions  in  an  externally-driven 
discharge  channel,  delivery  to  the  pellet  can  be  accomplished  with  the  time-varying  ion  energies 
appropriate  for  beam  bunching.  In  Sec.  V,  MHD  computations  illustrate  that  an  upper  limit  of 
about  1  MA  of  beam  current  exists  because  of  plasma  channel  expansion.  The  expansion 
reduces  magnetic  confinement  of  the  transporting  beam  causing  radial  beam  expansion  and  loss. 
Calculation  of  growth  rates  for  electrostatic  instabilities  in  the  beam-plasma  system  in  Section 
VI  F  indicate  that  this  class  of  instabilities  does  not  inhibit  transport  for  beams  satisfying  the 
MHD  constraints. 

Preliminary  light-ion  transport  experiments  are  discussed  in  Sec.  VI  G.  The  discharge 
channel  was  established  along  a  50  cm  thin  wire  embedded  in  a  reduced  air-pressure  back¬ 
ground.  Although  transport  was  demonstrated,  the  experiment  was  too  crude  to  determine 
transport  efficiencies.  Additionally,  the  transported  beam  current  was  too  low  to  assess 

difficulties  associated  with  MHD  response  of  the  plasma  channel  or  instabilities. 

Plasma  effects  encountered  by  focusing  heavy-ion  beams  traversing  reactor-chamber  back¬ 
grounds  are  summarized  in  Secs.  VI  H-J.  Calculations  of  the  degree  of  magnetic  neutralization 
of  the  beam  are  treated  in  the  first  two  entries.  These  calculations  are  important  because  even 
1%  non-neutralization  may  be  sufficient  to  disrupt  beam  focusing  on  target  for  highly-stripped 
heavy-ion  beams.  Regimes  for  electrostatic  microinstability  of  heavy  ion  beams  are  considered 
in  Sec.  VI  j.  As  with  light-ions,  stability  was  found  provided  the  beam  parameters  are  con¬ 
sistent  with  currently-conceived  reactor-level  beams. 
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II.  PARAVECTOR-POTENTIAL  DIODE  RESEARCH 

A.  INTRODUCTION  AND  THEORY 

Pinched  electron  flows  in  diodes  have  been  used  as  a  means  of  concentrating  relativistic 
electron  beam  kinetic  energy  onto  small  areas.  Inertial  confinement  fusion  with  relativistic 
electron  beams  has  been  a  prime  motivation  for  experimental  and  theoretical  efforts  directed  at 
maximizing  the  electron  power  concentration.  During  recent  studies,  it  was  found  that  hollow 
cathodes  generate  a  collapsing  hollow  ring1  that  pinches  because  of  an  ion-induced  pinch.2  3 
Prior  to  these  discoveries,  it  was  speculated  that  the  electron  flow  pinches  in  vacuum4  5  due  to 
the  self-magnetic  field  generated  by  the  electron  flow  current  when  it  exceeded  a  critical  current 
for  pinching.6  It  has  been  shown  that  the  relativistic  effects7  of  unneutralized  electron  flows 
under  self-magnetic  fields  do  not  generate  a  tightly  pinched  electron  flow,  but  rather  a  weakly 
pinched  one.  The  reason  for  this  is  that 

|E|sELF  >  l®lsELF  0) 

and  electrons  do  not  E  x  B  drift  under  such  conditions.  The  models4  5  for  vacuum  pinching, 

however,  can  be  made  self-consistent  by  adding  a  mysterious  current  flowing  on  the  diode  axis. 
This  current  adds  to  the  total  magnetic  field  reversing  inequality  (1).  The  resulting  electron 
flow  is  described  by  a  force-free  flow  which  has  been  coined4  5  Para-Potential  Flow  (PPF).  One 
may  then  summarize  the  ion-induced  pinch  and  the  PPF  described  above  for  generating 
pinched  electron  flow  as  two  ways  of  reversing  inequality  (1).  The  ion-induced  pinch  intro¬ 
duces  ions  that  reduce  the  electric  field  while  the  PPF  scheme  introduced  a  bias  current  that 
increases  the  magnetic  field.  The  latter  shall  be  called  a  bias-current  pinch  since  the  PPF  is 
only  one  of  a  family  of  electron  flows  that  ExB  drift  towards  the  diode  axis. 
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There  is  a  definite  advantage  of  bias-current  pinches  over  ion-induced  pinches  for  low 
impedance  (<  1  ft)  generators.  The  ion  current  that  flows  in  low  impedance  pinched-beam 
diodes  restricts2  8  the  electron  power  to  a  finite  level  (— 10I2W)  and  is  thus  an  unacceptable 
power  drain  for  electron-beam  fusion  (but  is  an  interesting  source  of  power  for  ion  beams2,9). 
Historically,  there  were  many  attempts  to  generate  bias-current  pinches.  These  included  the 
use  of  (1)  tapered  hollow  cathodes10,  (2)  a  plasma  column  on  the  diode  axis,"  (3)  an 
externally-driven  rod  on  the  diode  axis,12  (4)  two  concentric  cathodes,13  and  recently  (5)  multi¬ 
ple  cathodes.14 

Here,  we  report  on  a  new  technique  that  is  a  combination  of  a  bias-current  pinch  and  an 
ion-induced  pinch  (Fig.  1).  This  technique  allows  one  either  to  minimize  the  ion  current  or  to 
concentrate  it  near  the  axis.  In  addition,  electron  injection  into  the  diode  is  radial  rather  than 
the  axial  emission  of  electrons  from  cathodes  used  to  date.  The  radial  injection  was  attempted 
as  a  smooth  transition  for  electrons  from  the  cathode  to  the  equipotential  surfaces  (see  the 
electron  trajectories  shown  in  Figure  1).  Since  electrons  must  cross  equipotentials  starting  from 
the  cathode  (where  the  potential  <t>  =  0)  on  their  way  to  motion  along  the  equipotential  sur¬ 
faces  near  the  anode,  they  can  not  be  characterized  by  a  PPF  theory.  A  more  general  theory 
which  is  restricted  only  to  cold  electron  flow1  shows  that  the  mechanical  momentum  of  the 
electron  fluid  is  related  to  the  magnetic  field  by 

V  *  P  =  -  —  B  (2) 

c 

Since  B»VxA,  and  A  is  free  to  be  defined  with  the  addition  of  any  V*  one  may  realize 

P  =  —  —  A  (3) 

c 

by  choosing  V*  —  0.  Thus,  the  electron  mechanical  momentum  vector  lies  along  the  vector 
potential.  Hence  the  Paravector-potential  (PVP)  diode. 
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The  actual  solutions  of  the  electron  flow  for  the  geometry  of  Fig.  1  near  the  inner  edge  of 
the  outer  cathode  and  the  anode  edge  are  rather  complicated.  The  exact  solutions  at  smaller 
radii  will  depend  on  an  analytic  continuation  of  the  flow  pattern  from  the  regions  of  the  corners 
into  the  smooth  region  of  the  flat  anode.  A  family  of  solutions  which  are  exact  solutions  of  the 
generalized  equation7  is  given  by  a  treatment  similar  to  that  of  the  reference. 


The  solutions  are  based  on  the  existence  of  a  region  of  zero  electric  field  ((a)  in  Fig.  1), 
an  electron  sheath  of  thickness  d\  (region  (b)),  and  a  vacuum  region  of  thickness  d2  (region 
(c)).  The  magnetic  field  in  region  (a)  is  derived  from  the  bias  current  of  the  center  cathode, 
and  ions  flowing  to  it  from  the  anode.  Electrons  from  the  outer  cathode  flow  without  being  lost 
on  the  anode  at  radii  between  two  cathodes.  Using  the  forms  of  the  electric  and  magnetic  fields 
of  the  electron  flow  as  derived  by  Goldstein,  et  al.,7  the  thickness  of  the  electron  flow  du  and 
the  thickness  of  the  flow-anode  vacuum  gap  d2  are  found  as  functions  of  radius. 


here  r  is  the  radius,  /„  —  -  -  17  kA,  lc  is  the  center  current,  y„ 

e 

diode  voltage,  and  the  functions  Fand  G  are  given  by 


(4) 

(5) 


1  + 


eV, 

mec 2 


V0  the 


F(x)  -  1/2/n  lx  +  (x2  -  1)1/J] 
x+(x2-  1)'/J 


Gix )  -  1/2 


(6) 

(7) 


x2  +  x  (x2  —  1)1/J-  1 

Over  the  range  2  <  x  <  3,  Fix)  and  G ix)  are  approximately  0.8  and  0.25,  respectively.  It  is 
important  to  note  that  for  a  given  diode  voltage,  the  gaps  are  totally  defined  by  the  ratio  -j- 

*C 


The  anode-cathode  spacing  in  region  (a)  does  not  enter  into  the  result. 
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Since  d2  must  be  positive  (or  else  the  electrons  will  strike  the  anode),  lc  must  be  greater 

than  —  /,.  As  d2  is  the  distance  between  the  flow  and  the  anode,  it  provides  a  measure  of  the 

y0 

axial  anode  perturbations  that  the  diode  can  tolerate.  For  example,  for  diode  parameters  of  1 
MV,  I  MA,  and  y-  =  2.4,  the  quantity  d2  =  ,006r  while  d,  is  a  factor  of  5  greater.  Thus,  at  a 

*V 

radius  of  4  cm,  one  can  tolerate  anode  imperfections  on  the  order  of  0.02  cm. 

If  only  electron  flow  existed,  then  the  space  charge  effect  of  the  flow  from  the  outer 
cathode  as  it  came  opposite  the  inner  cathode  would  have  decreased  the  electron  flow  from  the 
inner  cathode,  thus  decreasing  the  bias  current,  resulting  in  the  loss  of  the  electrons  on  to  the 
anode  opposite  the  inner  cathode.  Taking  into  account  that  an  anode  plasma  is  being  produced 
opposite  the  inner  cathode  by  the  electron  beam  heating  of  the  anode  surface,  it  is  seen  that  ion 
flow  may  provide  space-charge  neutralization  negating  the  effect  of  the  radial  electron  flow  and 
replacing  the  partially-suppressed  electron  current.  The  sum  of  the  electron  and  ion  currents  in 
the  center  cathode  represents  the  needed  bias  current.  The  presence  of  the  ion  flow  allows  an 
ion-induced  pinch  to  continue  the  pinching  of  the  electron  flow  to  radii  smaller  than  the  inner 
cathode.  Numerical  simulation15  16  has  shown  that  in  this  region,  the  drifts  of  electrons  (such 
as  the  Vfi  drift)  are  a  prime  cause  for  finite  pinch  radius.  The  use  of  smaller-hole  cathodes  is 
known17  to  generate  smaller  pinches  and  was  also  tried  in  experiments  reported  on  here. 

In  the  numerical  simulation  studies,  an  extension  of  the  r  -  z  semistatic,  particle-in-cell 
code  was  used  to  find  the  self-consistent,  numerical  solutions  of  electron  and  ion  flows.  In  Fig. 
2,  the  diode  geometry  and  vacuum  equipotentials  are  shown  for  diode  parameters  close  to  those 
used  in  the  experiments.  When  self-consistent  particle  flows  are  generated  in  this  diode,  the 
equipotential  field  lines  (Fig.  3)  are  drastically  altered  and  shifted  toward  the  anode  in  close 
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agreement  with  analytic  'heory.  The  only  constraint  imposed  was  that  ions  on  the  anode  were 
generated  at  radii  smaller  than  that  of  the  inner  cathode. 

The  basic  concept  of  the  PVP  diode  led  to  the  design  of  a  series  of  diodes  used  in  con¬ 
junction  with  the  GAMBLE  I  and  GAMBLE  II  generators  operated  in  positive  polarity.  The 
diode  performance  was  studied  over  wide  variations  of  voltage  and  total-to-center  current  ratio 
/,//,.  In  agreement  with  theory,  it  was  found  that  given  a  total-to-inner  current  ratio  less  than 
y0,  the  radially-injected  electrons  will  flow  inward  without  striking  the  anode  at  large  radii  pro¬ 
viding  there  are  no  gaps  or  projections  to  perturb  the  equipotential  surfaces.  If  the  needed  bias 
current  is  not  provided,  then  a  self-correcting  phenomenon  was  observed  in  which  excess 
radially-injected  electrons  are  lost  immediately  to  the  anode  edge.  The  rest  flow  in  without 
striking  the  anode  at  large  radii  resulting  in  a  current  of  y„lt  near  the  anode  axis.  The  other 
predicted  phenomenon  observed  was  the  interaction  of  the  radial  electron  flow  with  the  bias 

flow.  Instead  of  continually  rising  during  the  voltage  pulse,  the  bias  current  clamped  at  the 
value  it  had  at  the  time  the  radial  flow  started. 

B.  EXPERIMENTAL  APPARATUS 

The  PVP  diode  design  has  worked  over  the  course  of  more  than  80  shots.  A  scale  draw¬ 
ing  of  the  latest  geometry  is  shown  in  Fig.  4.  The  cathodes  (shown  solid)  were  made  from 
brass  and  mounted  on  a  36  cm  diameter,  2.5  cm  thick  aluminum  insert  (shown  cross-hatched) 
which  fit  the  doors  of  both  of  the  GAMBLE  relativistic  beam  generators  when  operated  in  posi¬ 
tive  polarity  The  door  insert  has  a  4.4  cm  diameter  hole  on  axis  (originally  2.5  cm)  and  a  1.0 
cm  slit  22.9  cm  long  cut  through  it  (originally  19.7  cm)  for  viewing  a  diameter  of  the  anode. 
The  taper  of  the  inner  cathode  was  varied  from  6°  to  30°  with  a  variety  of  inner  and  outer 
diameters.  The  most  recent  inner  cathode  was  the  8. 2/1. 3  cm  (O.D./I.D  ),  12°  tapered  hollow 
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Figure  4  —  Scale  drawing  of  latest  PVP  diode 
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cathode  shown  in  Fig.  4.  The  outer  cathode  was  a  22.9  cm  I  D.,  2S.4  O.D.  brass  ring.  Shots 
using  annular  brass  inserts  to  close  the  radial  gap  or  tailor  the  emission  surface  have  been 
taken.  The  anodes  were  originally  flat  aluminum  disks  (Fig.  Sa)  20  cm  in  diameter  which  were 
held  onto  the  center  conductor  of  GAMBLE  I  or  GAMBLE  II  by  6  insulating  bolts  at  a  large 
radius.  On  later  shots,  the  anode  diameter  was  increased  to  close  the  radial  gap.  Spacers 
allowed  the  anodes  to  be  positioned  so  as  to  be  in  the  plane  of  the  outer  cathode.  Later,  in  an 
attempt  to  reduce  the  edge  damage,  the  outer  edge  was  beveled  at  45°  (Fig.  Sb)  to  give  the 
non-radially  emitted  electrons  more  distance  in  which  to  turn.  Only  in  the  last  series  of  shots 
was  the  outer  edge  geometry  changed  to  that  shown  in  Figs.  4  and  5c  —  the  anode  edge  being 
relatively  sharp  and  the  anode  beveled  at  negative  45°.  In  this  configuration,  the  axial  electric 
field  was  reduced  with  the  result  that  most  of  the  current  flowed  radially.  The  anode  was  a  20.3 
cm  diameter,  stainless  steel  outer  structure  (shown  with  fine  hatching)  which  held  an  aluminum 
insert  (shown  with  large  hatching).  The  space  behind  the  insert  was  present  to  prevent  shock 
damage  to  the  stainless  holder  and  the  hub  of  the  generator.  The  stainless  holder  was  mounted 
with  a  single  1/4"  —  20  screw  into  the  generator.  The  entire  anode  surface  was  covered  with 
either  2  ^m  aluminized  KIMFOL  a  polycarbonate  foil  or  6  (irn  aluminized  mylar,  aluminum 
side  towards  the  cathode.  The  covering  provided  the  anode  with  a  smooth,  continuous 
potential  surface. 

The  main  diagnostics  were  a  diode  voltage  monitor,  B  loops  monitoring  the  current 
flowing  from  the  center  cathode  and  from  the  total  cathode  assembly,  anode  witness-place  dam¬ 
age,  and  x-ray  production  monitored  by  a  time-integrated  x-ray  pinhole  camera  and  calibrated 
p-i-n  diodes. 
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C.  OBSERVATIONS 

Series  I  -  GAMBLE  I  Shots  #  5828-5836 

This  series  of  shots  used  20.0  cm  diameter  flat  aluminum  anodes  mounted  with  six  bolts 
at  a  diameter  of  18.4  cm.  The  bolts  were  counterbored  so  as  not  to  protrude  out  of  the  plane 
of  the  anode.  An  insert  into  the  outer  cathode  made  the  radial  gap  5  mm  (see  Fig.  5e).  The 
anode  was  varied  from  being  inserted  0.5  mm  into  the  plane  of  the  outer  cathode  to  being  with¬ 
drawn  0.5  mm  from  it.  The  inner  cathode  was  an  8. 4/3. 9  cm  (O.D/I.D.),  6°  tapered  hollow 
cathode  and  the  anode-cathode  gap  was  5.2  mm  at  the  inner  diameter. 

The  voltage  on  these  shots  was  low  (<  400  kV)  and  the  bias  current  criterion  was  not 
met.  Only  100  kA  out  of  the  235  kA  total  current  flowed  in  the  center  cathode  which,  accord¬ 
ing  to  the  theory  of  PVP  flow,  could  bring  in  80  kA  from  the  outer  cathode  leaving  55  kA  to 
strike  the  anode  rim.  Anode  rims  showed  existence  of  radial  flow  through  their  extensive  dam¬ 
age.  The  inner  edges  of  the  screw  holes  showed  evidence  that  the  beam  had  struck  them  on  its 
way  in  lending  support  for  the  existence  of  radial  flow. 

The  other  conclusion  from  this  series  of  experiments  is  that  the  very  large  area  anode 
(314  cm2)  could  not  turn  on  an  anode  plasma  during  the  50  ns  pulse  time  so  that  an  ion 
induced  pinch  could  not  take  place.  This  is  in  agreement  with  the  axial  velocity3  on  aluminum 
(~  2  cnt2/ns/100kA)  which  brought  the  collapsing  hollow  ring  from  a  radius  of  10  cm  to  a 
radius  of  9  cm. 


Series  II  —  GAMBLE  I  Shots  #  5837  and  5845 
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For  these  shots,  the  same  aluminum  anodes  were  used  as  well  as  the  same  inner  cathode. 
However,  the  outer  cathode  had  finger  stock  mounted  on  its  inner  edge  instead  of  an  insert 

(see  Fig.  5f) .  The  inner  cathode  provided  the  required  bias  current  and  was  nearly  equal  to  — 

To 

/,  for  almost  the  total  pulse  time.  The  anodes  showed  little  edge  damage  and  had  a  rear-surface 
spall  of  1.1  cm  diameter  on  axis.  The  x-ray  pinhole  photograph  for  shot  #5845  shows  a  7  mm 
(FWHM)  pinch.  Figure  6  shows  the  front  (left)  and  back  of  the  anode  for  shot  #5845.  Note 
the  slight  damage  at  the  inner  edge  of  the  screw  holes.  Figures  7  and  8  show  the  diode  param¬ 
eters  for  the  same  shot.  It  is  significant  to  note  that  the  outer  cathode  did  not  turn  on  for  10- 
20  ns  after  the  start  of  the  current  pulse.  When  the  outer  cathode  did  turn  on,  the  inner 

cathode  current  was  clamped  at  its  then  current  value  remained  nearly  equal  to  —  /,  for  the 

Vo 

rest  of  the  pulse.  The  value  at  which  the  inner  current  clamped  was  a  factor  of  two  below  its 
operating  value  when  used  without  the  outer  cathode.  This  agrees  with  our  theoretical  expecta¬ 
tions  of  the  effect  of  the  radial  electron  flow  on  the  behavior  of  the  inner  cathode. 

Series  III  -  GAMBLE  I  Shots  #5847-5851 

These  shots  were  identical  to  the  shots  in  Series  II  with  the  addition  of  a  TRW  streak 
camera  looking  at  the  ions  emitted  from  the  anode.  The  ions  struck  the  silvered  side  of  an 
optically-thin  slab  of  Pilot-B  scintillate r  covering  a  diameter.  The  scintillator  was  viewed  with 
the  streak  camera.  The  streak  camera  photographs  showed  that  there  were  no  ions  emitted  at 
large  radii  (compared  to  the  outer  radius  of  the  inner  cathode)  although  there  were  some  ions 
coming  from  just  outside  the  inner  cathode. 
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Figure  6  —  Aluminum  anode  for  shot  #5845 
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Series  IV  -  GAMBLE  I  Shots  #  5861-5868 

These  shots  had  a  45°  beveled  outer  cathode  insert  (Fig.  5g)  and  a  20.3  cm  or  20.6  cm 
diameter,  45°  beveled  anode  (Fig.  5c).  The  reduced  radial  gap  meant  that  the  voltage  was  again 
low  and  the  required  bias  current  not  supplied.  This  geometry  did  not  prove  successful. 

Series  V  -  GAMBLE  II  Shots  #1731-1735 

For  the  first  series  of  shots  on  the  higher  power  GAMBLE  II  generator,  20.6  cm  diame¬ 
ter,  45°  beveled  anodes  (Fig.  5c)  and  an  outer  cathode  with  finger  stock  mounted  on  its  inner 

edge  (Fig.  5f)  were  used.  On  none  of  these  shots  was  the  —  /,  bias  current  supplied. 

y0 

Nevertheless,  shot  #1735  appeared  to  be  a  relative  success. 

Figures  9  and  10  show  the  diode  parameter  for  shot  1735.  The  diode  operated  at  750  kV, 
900  kA  with  a  flat,  0.50  impedance  during  the  plateau  of  the  voltage  pulse.  Two  things  should 
be  noticed.  1  he  bias  current  required  to  vacuum  pinch  (350  kA)  all  of  the  outer  current  was 
never  supplied  and  there  was  no  delay  time  between  the  start  of  the  inner  and  outer  currents. 
Even  so,  400  kA  out  of  a  outer  current  of  650  kA  should  have  vacuum  pinched. 

Figure  11  shows  the  front  (left)  and  back  of  the  anode  plate  for  shot  #1735.  Although 
there  was  a  great  deal  of  damage  to  the  anode  edge  due  to  the  250  kA  of  excess  current,  most 
of  the  anode  outside  of  the  inner  catnode  radius  is  relatively  undamaged.  Most  of  the  apparent 
damage  seen  in  Fig.  1 1  is  aluminum  melted  by  the  pinch  and  deposited  on  the  anode  surface. 
The  x-ray  pinhole  photograph  confirms  this  as  it  shows  relatively  little  x-radiation  coming  from 
the  region  between  the  anode  edge  and  the  inner  cathode  compared  to  radiation  from  the  anode 
edge.  This  provides  additional  support  for  our  theoretical  expectation  that  ynlc  will  pinch  with 
the  excess  outer  current  sloughing  off  to  the  anode  edge. 


24 


NRL  MEMORANDUM  REPORT  4237 


Figure  9  —  Aluminum  anode  for  shot  #1735 
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Series  VI  -  GAMBLE  II  Shots  #1745-1748 

These  shots  were  identical  to  the  last  series  except  that  the  outer  cathode  had  no  finger 
stock  insert  (Fig.  5d)  in  order  to  open  the  radial  gap,  and  the  inner  cathode  was  a  2.S  cm  I.D., 
2.9  cm  O.D.,  30°  tapered,  hollow  cathode.  The  sharp-edged  center  cathode  tended  to  either  not 
supply  very  much  current,  or  to  short  out.  However,  for  shot  #1748,  the  biasing  criteria  was 
satisfied  through  the  peak  of  the  voltage  pulse  (I  MV)  but  was  not  satisfied  during  the  voltage 
fall.  Again,  the  x-ray  photograph  shows  little  radiation  between  the  anode  edge  and  the  inner 

cathode. 

Series  VII  -  GAMBLE  I  Shots  #  6431-6447 

This  last  series  of  shots  used  the  geometry  shown  in  Fig.  4,  the  design  of  which  was  based 
on  considerations  described  earlier.  The  electrical  characteristics  and  x-ray  photographs  showed 
that  this  new  geometry  almost  always  met  the  biasing  criterion  and  that  there  were  few  or  no 
x-rays  from  the  anode  edge.  Also,  this  diode  operated  with  a  shot-to-shot  reproducibility  of 
10%  in  current  and  voltage.  Figures  12  and  13  are  the  averaged  diode  behavior  for  the  five 
shots  6442-6446.  As  was  typical  of  the  earlier  successful  GAMBLE  I  shots,  there  is  a  10-20  ns 
time  delay  between  the  start  of  the  total  current  and  the  start  of  the  outer  current. 

In  order  to  study  the  ion  current,  the  anode  was  modified  by  drilling  a  1.3  cm  diameter 
hole  on  axis  and  mounting  a  2  cm  diameter  disk  of  127  *im  polyethylene  over  the  hole.  The 
anode  was  recessed  slightly  over  the  central  2  cm  so  that  the  polyethylene  disk  would  not  pro¬ 
trude.  The  mylar  anode  covering  was  cut  away  to  expose  this  disk.  The  result  of  activation 
measurements  on  this  system  indicates  that  ion  current  densities  of  approximately  4  kA/cm 2 
flow  in  the  2.5  cm  diameter  hole  of  the  center  cathode.  These  preliminary  results  will  be 
further  investigated. 
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Figure  13  —  Average  Inner.  Outer,  Total,  and  Theoretical  Bias  Current  for  shots  #6442-6446 
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D.  CONCLUSIONS 

A  comparison  of  analytic  theory,  numerical  simulation  results,  and  experiments  indicates 
that  bias-current  pinching  via  paravector-potential  flow  may  indeed  be  realized.  Detailed 
geometries  of  diodes  specifically  designed  to  allow  a  smooth  transition  from  a  PVP  flow  at  large 
radii  to  an  ion  induced  pinch  flow  at  small  radii  were  studied.  With  careful  optimization,  it  may 
be  possible  to  create  very  efficient  electron  pinches  from  large  aspect  ratio  diodes. 

It  has  been  shown  that  the  PVP  diode  concept  allows  the  multiplication  of  the  total 
current  of  a  single,  tapered,  hollows  cathode  by  a  factor  of  y0  without  changing  the  geometry  of 
the  hollow  cathode.  In  addition  a  fraction  (l-l/y0)  of  this  current  may  be  electron  current 
independent  of  the  diode  aspect  ratio  or  impedance.  This  contrasts  with  single,  hollow  cathodes 

which  rely  on  ion  induced  pinching  and  contain  a  smaller  fraction  of  their  total  current  in  elec¬ 
tron  current.  The  same  diode  may,  however,  operate  with  a  bias  current  exceeding  —  /,  by 

y« 

allowing  efficient  ion  production  inside  the  center  cathode.  Preliminary  results  are  encouraging 
and  will  be  investigated  further. 
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III.  ELECTRON  TRANSPORT  EXPERIMENTS 

A.  INTRODUCTION 

The  propagation  of  relativistic  electron  beams  in  plasma  channels  has,  in  the  last 
decade,1-7  received  considerable  interest  due  to  its  applications  for  energy  transfer  and  pellet 
fusion.  A  rough  division  of  this  effect  may  be  made  by  characterizing  the  electron  beam  as  hot 
or  cold.  In  both  cases,  the  relativistic  electron  beam  propagates  in  a  plasma  having  a  self- 
consistent  azimuthal  magnetic  field  which  confines  the  electrons  radially.  The  necessary  mag¬ 
netic  field  for  hot  beam  propagation  (electrons  injected  with  a  large  mean  angle)  is  given  by 
Alfven’s  treatment  of  electron  orbits  which  yields  the  condition  that  the  net  current  inside  the 
beam  radius  must  be  in  the  same  direction  as  the  beam  current  and  nearly  equal  to  the  Alfven 
current  for  the  energy  of  the  relativistic  electrons. 

/,*  =  -  17  (y2  -  1)'/2  kA  (1) 

On  the  other  hand,  for  relatively  cold  beams  (mean  angle  of  injection  9  «  1),  the  same  treat¬ 
ment  indicates  that  the  current  needed  is  only 

/,„  -  (1  -  cos  9)  /„  a  y 92l„  (2) 

In  order  to  propagate  1  MeV  electrons  with  a  10°  mean  injection  angle,  one  needs  only 
about  1  kA  of  net  current  in  the  channel  compared  to  the  50  kA  of  net  current  needed  for  hot 

beam  propagation.  For  typical  beam  currents  of  a  few  hundred  kiloamperes  injected  into 
plasmas  at  densities  of  the  order  of  1018  cm-3,  the  return  current  induced  in  the  plasma  is 

nearly  equal  to  the  primary  beam  current,  leaving  a  net  current  of  the  order  of  1%  of  the  beam 
current.  The  remainder  of  the  return  current  nows  outside  of  the  beam  radius.  The  resulting 
net  current  of  a  few  kiloamperes  within  the  beam  radius  is  sufficient  to  propagate  relatively  cool 
beams  (as  discussed  above)  but  is  insufficient  for  the  propagation  of  a  hot  beam.  Other  means 
for  providing  the  needed  net  current  must  therefore  be  utilized  for  hot  beam  propagation.  His¬ 
torically,  these  included  z-pinch  plasmas5,  and  axial  discharges  via  exploding  wires.6  7 
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In  the  present  work,  the  physics  of  the  propagation  of  hot  beams  in  plasma  discharges 
driven  by  exploding  wires  is  described.  The  work  extends  previous  work6  7  into  a  regime  of 
higher  beam  currents  (up  to  300  kA  injected)  and  investigates  the  effect  of  channel  current  and 
diode  voltage  on  beam  propagation.  The  effects  of  scattering  and  axial  electric  fields  on  elec¬ 
tron  orbits  as  well  as  the  limitations  on  the  maximum  beam  current  and  beam  current  density 
that  can  be  propagated  is  investigated. 

B.  EXPERIMENT  A  L  APPAR  A  TUS  A  ND  OBSER  VA  TIONS 

The  experimental  results  of  the  present  work  were  obtained  on  the  GAMBLE  I  (800  kV, 
220  kA,  80  ns  FWHM)  and  GAMBLE  II  (1000  kV,  420  kA,  50  ns  FWHM)  electron  beam 
facilities.  A  schematic  of  the  experimental  apparatus  is  shown  in  Fig.  1.  An  84/39  mm 
(O.D./I.D.)  hollow  cathode  with  a  6°  taper  and  5.1  mm  anode-cahode  gap  (at  the  inner  diame¬ 
ter)  was  used.  The  anode  was  a  1.5  mm  thick  brass  plate  with  a  20  mm  diameter  hole  on  axis. 
The  hole  was  beveled  at  30°  on  the  cathode  side  and  sanded  to  provide  a  smooth  transition 
region  for  the  radial  electron  flow.  The  hole  was  covered  with  a  25  fi m  titanium  foil  which  was 
both  ground  plane  and  vacuum  window.  Immediately  behind  this  was  mounted  a  20  jum  alumi¬ 
num  foil  to  which  a  51  #xm  tungsten  wire  was  affixed.  The  wire  was  stretched  in  air  from  the 
center  of  the  foil  through  the  aluminum  target  plate,  located  30  to  64  cm  away,  to  a  56  /uf,  2.4 
jih  (total  circuit  inductance),  20  kV  capacitor  bank.  A  "squirrel  cage"  of  aluminum  rods  and 
rings  served  as  a  symmetric  current  return  for  the  channel  current.  This  bank  could  drive  a 
peak  current  of  50-80  kA  (occurring  approximately  15  *is  into  the  pulse)  in  an  air  plasma  sur¬ 
rounding  the  wire 

The  experimental  diagnostics  included  diode  voltage  and  current  monitors,  time  integrated 
hard  x-ray  pinhole  cameras  looking  at  the  anode,  plasma  channel,  and  target,  open  shutter 
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Figure  I  —  \  schematic  of  the  experimental  apparatus  showing  at  the  brass  cathode,  bt  the  brass  anode,  cl  the 
ntanium/aluminum  foil  assembly,  d)  the  tungsten  wire.  e>  the  aluminum  target.  *;id  ft  the  switch  and  capacitor  bank 
which  supplies  the  channel  current 
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visual  photography  of  the  channel,  time  resolved,  calibrated  x-ray  p-i-n  diode  detectors,  and 
self-integration- So  and  Rogowski  coils  monitoring  the  channel  current.  In  addition,  thermo¬ 
luminescent  detectors  were  used  to  obtain  the  peak  energy  of  the  transported  electrons.  Streak 

camera  measurements  of  the  plasma  channel  over  the  times  of  interest  allowed  one  to  deter¬ 
mine  the  channel  radius  at  the  time  of  beam  injection.  This  was  verified  by  the  time-integrated 
x-ray  photographs  of  the  channel.  The  channel’s  resistance  was  measured  by  the  shift  of  the 
frequency  of  the  ringing  channel  current  from  its  zero  resistance  value  -  measured  by  replac¬ 
ing  the  wire  with  a  1.5  cm  diameter  conducting  rod  -  and  by  the  current  decay  time.  Using  an 
average  channel  radius  ( rch  =  1.6  cm)  allowed  one  to  estimate  the  channel's  conductivity  to  be 
<r  ~  2  x  1014  s  ',  and  temperature  to  be  T  =  2  eV.  These  parameters  are  consistent  with 
theory.8 

On  GAMBLE  I,  the  channel  current  was  varied  from  about  y  /„  to  2  /„,  while  the  anode 

hole  (which  limited  the  maximum  injection  radius)  was  varied  from  less  than  to  greater  than 
the  channel  radius.  Propagation  was  successful  only  when  the  channel  and  beam  currents  had 
the  same  direction,  in  agreement  with  the  theoretical  prediction.  In  addition,  it  also  became 
obvious  that  a  crucial  element  of  the  experiment  was  the  diode-channel  interface.  It  was 
observed  that  any  asymmetrical  magnetic  field  al  the  target  location  caused  a  distortion  in  the 
electron  distribution  at  the  target.  The  inference  to  a  similar  effect  at  the  diode-channel  inter¬ 
face  is  made  even  more  important  as  any  current  feed  assymetry  at  the  diode-channel  interface 
could  result  in  stray  magnetic  fields  in  the  diode  which  could  affect  not  only  the  beam  injection, 
but  also  the  diode  physics. 

Table  1  documents  some  important  GAMBLE  1  shots.  Shown  are  the  diode  voltage  and 
current,  the  channel  current  and  radius  (determined  from  the  beam  injection  time  and  the 
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plasma  expansion  velocity— about  0.8  cm//is— and  confirmed  by  x-ray  and  visual  streak  photog¬ 
raphy),  the  radius  ra  which  enclosed  the  Alfven  current,  the  radius  of  front  surface  damage  rn , 
and  the  radius  of  the  back  surface  spall  rs. 


Table  1  —  Experimental  Parameter  Summary 


Shot  # 

F(kV) 

/(kA) 

/,h(kA) 

rth(cm) 

ra  (cm) 

rD  (cm) 

rs(cm) 

5693 

700 

230 

80 

1.6 

1.1 

1.1 

— 

5695 

700 

230 

77 

1.4 

5700 

800 

220 

68 

1.6 

1.2 

1.2 

5701 

700 

185 

58 

1.6 

1.3 

1.4 

5903 

510 

205 

42 

1.4 

1.2 

1.2 

5955 

540 

270 

38 

1.0 

_ 

0.9 

0.9 

0.6 

Figure  2  shows  the  diode  voltage  and  current  and  the  target  x-rays  for  Shot  5700  on 
GAMBLE  I.  Thermoluminescent  detectors  allowed  one  to  obtain  an  electron  kinetic  energy 
loss  in  the  anode  foil  and  plasma  channel  of  80  keV.  Using  this  energy  loss  and  the  diode 

voltage  gave  the  energy  of  the  electrons  striking  the  target.  The  collimated  p-i-n  detector 
looked  at  the  central  3.5  mm  radius  of  the  target.  From  the  electron  energy  and  the  p-i-n  sig¬ 
nal,  a  peak  central  current  of  22  kA  was  calculated.  Extrapolating  the  inferred  current  density 
out  to  the  spall  radius  implied  a  propagated  current  of  145  kA.  This  is  out  of  195  kA  of  elec¬ 
tron  current  in  the  diode  (220  kA  total  diode  current  -  25  kA  ion  current  computed  using  the 
Goldstein-Lee  formulag)  or  a  propagation  efficiency  of  74%. 

On  GAMBLE  II,  no  parameter  study  was  done.  The  goal  was  to  propagate  a  high  beam 
current  density.  On  four  shots,  calibrated  x-ray  p-i-n  detectors  indicate  that  more  than  200  kA 
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of  electrons  propagated  in  the  channel.  In  addition,  a  fifth  shot  is  believed  to  have  propagated 
more  than  200  kA  based  on  a  comparison  of  target  damage  with  known  shots. 

On  GAMBLE  II,  the  wire  channel  was  independently  viewed  with  another  p-i-n  detector. 
It  indicated  that  less  than  20  kA  of  electrons  were  lost  in  the  plasma  channel.  This,  however,  is 
an  upper  limit  as  it  assumed  that  the  primary  collisions  occurred  between  the  electrons  and  the 
low-Z  plasma  channel  rather  than  the  tungsten  wire. 

Figure  3  shows  the  diode  voltage  and  current  for  GAMBLE  II  Shot  1799.  In  addition,  the 
calibrated  p-i-n  signal  is  shown.  Its  peak  signal  corresponds  to  a  current  of  approximately  160 
kA.  Since  the  p-i-n  viewed  only  the  center  5.3  mm  radius,  this  indicated  an  average  current 
density  over  this  region  of  about  180  kA/cm2.  From  the  x-ray  pinhole  photograph,  the  area  of 
the  target  x-rays  (FWHM)  was  determined  to  be  approximately  1.6  cm3  and  the  total  pro¬ 
pagated  current  to  be  about  280  kA.  This  is  out  of  a  diode  electron  current  of  330  kA  (420  kA 
total  diode  current  —  90  kA  ion  current)  or  an  efficiency  of  approximately  85%.  Two  other 
shots  gave  propagation  efficiencies  over  50%,  one  of  which  had  the  p-i-n  looking  at  the  whole 
target  plate  and  giving  255  kA. 

Figures  4  and  5  show  the  x-ray  pinhole  photographs  for  the  shots  described  above.  At  the 
top  of  Fig.  4  is  a  photograph  taken  with  two  different  pinholes;  the  lower  providing  greater 
resolution.  (Since  the  wire  itself  appeared  as  a  line  source,  the  size  of  its  image  provided  a 
measure  of  the  resolution.)  On  the  left,  electrons  striking  the  edge  of  the  hole  in  the  brass 
anode  produced  the  annulus  surrounding  the  x-ray  image  of  the  pinched  beam  passing  through 
the  titanium  anode  foil.  In  the  top  image,  the  x-rays  produced  by  electrons  in  the  plasma  chan¬ 
nel  can  be  seen.  The  target  is  not  seen.  At  the  bottom  of  Fig.  4  is  a  photograph  taken  with 
another  pinhole  camera  which  shows  the  entire  system  from  left  to  right.  Since  the  target  was 
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Figure  4  —  X-ray  pinhole  photographs  for  GAMBLE  I  Shot  3700 
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II 


Figure  5  —  X-ray  pinhole  photograph  for  GAMBLE  II  Shot  1799 
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closer  to  the  camera  than  the  anode,  the  relative  magnifications  are  different.  On  the  left  is  the 
anode  with  the  inner  halo  corresponding  to  an  end-on  view  of  the  plasma  channel.  At  the 
center  of  the  anode  is  an  image  of  the  injected  pinched  beam  while  at  right  can  be  seen  the 
image  of  the  electrons  as  they  strike  the  target.  Figure  5  shows  the  propagation  of  electrons  for 
GAMBLE  II  shot  1799.  The  circular  ring  around  the  target  is  pinhole  shine-through  while  the 
crosshairs  used  for  alignment  are  seen  as  an  absence  of  exposure. 

C.  DISCUSSION 

The  problem  of  relativistic  electron  flow  in  plasmas  heated  before  beam  injection  has  been 
treated  at  length  in  the  literature.10  Most  of  the  treatments  were  concerned  with  the  plasma 
current  response  to  the  injected  beam  current.  Two  effects  were  neglected. 

First,  the  conducting  plasma  has  a  finite  radius  outside  of  which  there  exists  only  low  con¬ 
ductivity,  cool  air  (theoretical  models  assumed  high  conductivity  at  all  radii).  The  total  return 
current  is  forced  to  run  in  the  conducting  plasma  due  to  the  voltage  increase  on  the  target  for 
the  relativistic  electron  beam.  This  differs  from  present  theory  which  includes  inductive  effects 
only.  The  fact  that  the  return  current  is  forced  to  run  in  the  plasma  also  explains  why  the 
loop  placed  outside  the  plasma  (in  both  NRL  and  SLA  experiments)  showed  essentially  zero 
current  change  during  beam  transport.  There  may  be  about  0.1%  of  the  primary  beam  current 
going  through  the  capacitor  bank  and  the  return  current  rods  connected  to  the  anode  plane  of 
the  generator"  but  this  is  too  small  a  current  to  be  delected  with  the  loop  in  the  present 
experiments. 

Second,  the  hydrodynamics  of  the  plasma  has  not  received  proper  treatment.  Usually, 
magnetic  fields  are  viewed  as  confining  fields  perpendicular  to  their  direction  because  of  their 
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pressure  tensor.  In  the  case  of  beam  injection  into  a  conducting  plasma,  however,  the  magnetic 
field  acts  to  confine  the  beam  particles  but  at  the  same  time  the  plasma  is  pushed  outwards. 
The  larger  the  magnetic  field  the  greater  is  this  expansion  force.  The  reason  for  the  reversal  of 
the  role  of  the  magnetic  field  is  that  the  return  current  set  up  in  the  plasma  in  response  to  the 
injected  beam  flows  in  the  direction  opposite  to  the  beam  current.  This  plasma  current  jp  = 
-  jb  experiences  a  jp  xBt  force  that  causes  plasma  expansion.  For  constant  beam  current  den¬ 
sities  and  for  times  shorter  than  magnetic  diffusion  times,  the  force  causes  constant  acceleration 
with  a  radial  expansion  proportional  to  the  time  squared  up  to  the  beam  pulse  duration. 


__  Jp  *  r2  __  IQ-12  WehT 


where  the  beam  and  channel  currents  Ub,  lch)  are  given  in  kA,  beam  pulse  time  in  10~8  s, 
plasma  radius  is  assumed  equal  to  beam  radius  and  given  in  cm,  and  the  plasma  density  p  is 
given  in  g/cm3.  The  expansion  is  less  that  10~2  cm  for  the  present  and  past6  7  experiments 
because  of  large  radii  (r  —  1  cm)  and  low  current  (~  105  A).  If  however,  106  A  would  have 


been  attempted  with  r  =  10-1  cm,  the  plasma  would  have  expanded  more  than  its  initial  radius 
during  50  ns  and  the  current  density  would  have  dropped  to  the  105  A/cm2  level. 


It  is  thus  concluded  that  the  present  experiments  in  open  air  were  performed  at  high  den¬ 
sity  with  little  hydrodynamic  effect  during  the  beam  propagation  while  Physics  International 
experiments5  with  tapered  Z-pinches  and  low  densities  ip  ~  10~7  g/cm3)  were  strongly  affected 
by  hydrodynamic  expansion. 

Some  limitations  on  current  propagation  are  now  considered.  The  simplest  limit  is  that 
due  to  the  plasma  conductivity.  An  electric  field  £  -  jb/<rp  is  needed  to  drive  jp  =  -  jb.  An 
upper  limit  on  this  field  for  one  Gamble  II  experiment  is  16  kV/m.  The  question  of  magnetic 
field  diffusion  and  its  effect  on  electron  beam  trajectories  is  now  considered. 
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Because  of  finite  channel  conductivity,  the  magnetic  field  profile  will  change.  The  time 
scale  for  change  over  a  length  r  in  the  radial  direction  is  given  by 

4ir<rr2 

'  =  -71- 

For  a  length  scale  r  —  0.5  cm  and  plasma  conductivity  of  cr  =  2  x  1014  s~\  t  —  750  ns. 
The  plasma  temperature  of  2  eV  is  obtained  from  heating  by  the  discharge  current  during  20 
fis.  Little  additional  heating  is  expected  for  H^A/cm2  beam  currents  and  the  plasma  densities  in 
the  experiment.  During  a  100  ns  pulse,  the  net  current  will  then  increase  by  A /  =  0.15  /*,  or 
30  kA  for  lb  =  200  kA.  The  increased  magnetic  field  due  to  the  net  current  drastically  affects 
the  electron  orbits.  If  the  beam  radius  at  injection  is  nearly  equal  to  the  plasma  channel  radius, 
the  relativistic  electron  injected  at  large  angles  will  be  returned  to  the  diode  by  this  magnetic 
field  and  the  efficiency  of  propagation  will  be  reduced.  If,  however,  the  beam  radius  is  initially 
smaller  than  the  channel  radius,  then  the  electrons  will  expand  radially  to  a  radius  smaller  than 
the  channel  radius.  As  the  net  current  increases,  this  radius  will  decrease  and  the  beam  will 
propagate  at  a  higher  current  density.  No  reduction  in  efficiency  will  then  be  observed.  This 
effect  should  be  more  pronounced  at  higher  total  beam  currents.  The  experiments  on  Gamble 
I,  Hydra  and  Gamble  II  are  characterized  by  beam  currents  of  100,  200,  and  250  kA  respec¬ 
tively.  In  view  of  the  above  discussion,  one  should  have  observed  the  tightest  pinches  pro¬ 
pagated  on  Gamble  II.  This  was  indeed  observed. 

The  effect  of  time  variation  of  diode  voltage  is  important  since  a  channel  current  nearly 
equal  to  /„  is  needed  for  transport.  When  diode  voltage  is  reduced  the  needed  channel  current 
is  accordingly  reduced.  The  combination  of  initial  channel  current  and  increase  by  magnetic 
diffusion  brings  the  total  current  above  /„  as  diode  voltage  falls.  This  is  one  of  the  reasons  why 
Gamble  I  or  Hydra  (which  had  nearly  constant  voltage  pulses)  show  almost  constant  efficiency 
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of  transport,  while  fast-changing  impedance  experiments  on  Gamble  11  showed  a  drastic  reduc¬ 
tion  of  beam  transport  efficiency  during  the  SO  ns  pulse. 

Based  on  the  preceding  discussions,  limits  on  the  total  current  that  can  be  transported 
using  the  present  technique  can  now  be  determined.  First,  the  conductivity  defines  an  electric 
field  which  cannot  exceed  the  diode  voltage  divided  by  the  channel  length.  For  V  —  1  MV, 
a  —  2  x  1014  s'1  /  «  1  m,  the  limiting  current  density  is  2  x  106  A/cm2.  For  a  distance  of  10 
meters,  only  2  x  105  A/cm2  can  be  transported  in  these  channels.  At  the  higher  current  densi¬ 
ties  estimates  will  be  modified  by  the  higher  temperatures  and  conductivities  associated  with 
beam  and  return-current  heating. 

Second,  the  plasma  channel  expansion  limits  the  tightness  of  the  transported  beam.  If  the 
2  MA/cm2  beam  of  1  MV  electrons  discussed  above  is  injected  into  a  5  mm  radius  air  channel 
carrying  the  Alfven  current,  the  channel  radius  will  double  in  50  ns.  Even  taking  magnetic 
diffusion  into  account,  the  Alfven  radius  will  still  be  1.2  times  the  original  channel  radius. 

It  has  been  shown  experimentally  that  high-power-density  electron  beams  may  be 
efficiently  transported  using  current-carrying  air-plasma  channels.  It  has  been  shown  that  cer¬ 
tain  restrictions  apply  to  their  use,  the  most  restrictive  being  the  limitation  on  the  current  den¬ 
sity.  This  restriction  may  be  overcome  to  a  certain  extent  by  increasing  the  channel  tempera¬ 
ture.  This  might  be  accomplished  by  using  a  lighter  background  gas.  The  self-consistent  prob¬ 
lem  of  plasma-channel  expansion  occurring  simultaneously  with  the  magnetic  diffusion  and  its 
ultimate  effect  on  the  transported  beam  and  the  transport  efficiency  must  still  be  studied. 
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IV.  REFLEXIVE  FOCUSSING 

Purely-ballistic  focussing  is  not  appropriate  in  a  light-ion  driven  fusion  reactor  scenario 
where  large  focal  distances  are  required  This  is  because  of  limits  imposed  by  transverse  ion 
motion  brought  about  by  finite  anode  temperature,  scattering,  or  fluctuating  diode  fields.  If 
these  effects  lead  to  a  transverse-thermal-velocity  spread  in  the  beam  characterized  by  V,h,  then 
the  minimum  size  to  which  the  beam  may  be  focussed  can  be  approximated  by 

a  =  LAff  (1) 

where  L  is  the  focal  distance  from  the  diode  to  the  target  and  Aft  is  the  beam  divergences 

=  VJ  f  „.  Here,  V„  is  the  axial  ion  velocity.  Now,  if  a  —  1  cm  and  L  ~  103  cm  as  required 

for  a  reactor,  allowable  values  of  Aft  are  unrealistically  small  by  one  to  two  orders  of  magnitude. 

(Based  on  current  knowledge.  Aft  in  the  range  of  a  few  degrees  appears  reasonable.)  One  must 

then  seek  a  magnetically-assisted  focussing  technique  to  provide  the  necessary  standoff. 

Magnetic-focussing  techniques  have  been  discussed  in  the  literature.  However,  neither 
technique  is  appropriate  for  the  delivery  of  Af4/cm2  ion  beams  to  distant  targets.  Both  involve 
the  use  of  large  applied  magnetic  beams  to  distant  targets.  Both  involve  the  use  of  large  applied 
magnetic  fields  in  the  diode  and  are  therefore  appropriate  to  a  reflex-triode  mode  of  ion  beam 
production.  The  low  ion-current-density  capability  of  this  technique  coupled  with  an  upper 
limit  of  an  order-of-magnitude  increase  in  ion-current  density  associated  with  the  two  focussing 
schemes  leads  to  to  deliverable  ion-current-densities  which  are  too  low  to  meet  pellet  require¬ 
ments.  Although  self-pinch  diodes  can  produce  the  required  current  density,  the  applied  field 
must  be  excluded  from  the  diode  for  the  ion  production  mechanism  to  operate.  In  that  case, 
the  large  angular  velocity  induced  as  the  emerging  beam  enters  the  field-filled  region  then 
prevents  focussing.  Additionally,  both  focussing  techniques  employ  curved  magnetic-field  lines 
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which  induce  currents  in  the  propagating  charge-neutral  beam.  In  order  that  these  currents  not 
disrupt  the  focus,  the  energy-density  of  the  beam  cannot  exceed  that  of  the  applied  fluid;  a  con¬ 
dition  which  cannot  be  satisfied  in  the  vicinity  of  the  focus. 

For  the  problem  at  hand,  it  may  be  desirable  to  carry  out  focussing  with  standoff  by 
employing  weak  guide  fields  with  straight  lines  of  force.  This  can  be  accomplished  by  introduc¬ 
ing  the  beam  at  a  small  angle  to  the  field  and  allowing  the  large  gyro-radius  motion  of  the  ions 
to  bring  the  beam  to  focus  meters  downstream.  Since  the  ions  must  move  across  the  magnetic 
field,  charge  and  current  neutralization  cannot  be  provided  by  electrons  extracted  from  the 
diode.  A  background  plasma  or  gas  fill  is  therefore  required.  The  reflexive-focussing  technique 
just  described  is  shown  schematically  in  Fig.  1.  The  anode  is  curved  so  that  the  beam  is  made 
to  diverge.  The  symmetry  of  the  single-particle  motion  about  the  midplane  causes  the  real 
focus  on  the  right-hand-side  of  the  figure  to  reflect  the  virtual  focus  to  the  left  to  the  diode. 
Thus,  although  transported  over  a  long  distance  T,  the  final  focus  is  characterized  by  the 
shorter  length  For  this  reason,  thermal  beam  spreading  does  not  severely  limit  Focusabilily 
since  L  is  replaced  by  r,  in  Eq.  (1).  The  remainder  of  this  section  is  devoted  to  a  kinetic- 
theory  analysis  of  this  concept. 


It  is  assumed  that  complete  charge  and  current  neutralization  exists  so  that  the  ions  move 
only  in  the  applied  solenoidal  field.  The  beam-ion  distribution  function  therefore  satisfies 
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where  Jl  =  cB/m  is  the  ion  gyrofrequency.  Forms  for  /  which  depend  only  on  the  single¬ 
particle  constants  of  the  motion  satisfy  this  equation,  i.e. 


/  =  /(£.  Kp.U) 


(3) 
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where 

l  -  r.;  E  —  r:  -  ry;  p  =  n-„  +  V;  1)  r:: 

L  =  11  rv.  cos  £  -  (ly  +  v,;  +  rilvw)(sin<;  —  sin  d>)  +  ' :  11  V  sind) 


14) 


and  ;  =  ll  r/ v  +■  <t>.  The  quantity  <b  is  a  constant. 

The  desired  solution  must  satisfy  a  number  of  requirements.  The  distribution  function 
must  reflect  the  radial  divergence  of  the  beam  at  injection.  The  beam  must  be  allowed  to  haw 
rigid  rotation  within  certain  limits  at  injection  since  all.  or  any  fraction  of.  of  the  solenoidal  field 
may  be  allowed  to  penetrate  the  diode  region.  Since  the  beam  is  born  with  zero  angular 
momentum,  the  beam  emerging  inti'  the  drift  section  will  have  art  angular  velocity  (gi\en  by 
conservation  of  />)  of 


i..i-i  M!  !!>'  w 

where  11,,  is  the  ion  g>  rofrequeiuv  awa.iusl  with  the  'olenmd.il  field  in  the  diode  lhus. 
~U  2  ^  in  0  tor  positive  11  Inc  he  am  i.tdu  s  must  he  finite  (charactetized  hv  tadius  R  at 
injection)  and  a  transverse  ihcnoal  ve  im;i  spiead  is  desited  1  malls .  the  distribution  (unction 
must  be  analytically  tractable  \n  appropriate  form  lot  t  is 
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where  //(\)  =  1  for  \  sj  0  and  zero  otherwise  The  quantity  (  therefore  specifies  the  beam 
radius  at  injection  into  the  drift  space  lr  =  ()).  and  n  is  the  axial  beam-ion  density  there 


Although  satisfying  all  of  the  above  requirements,  1  q  (  ')  relaxes  exact  energy  conserva¬ 
tion  tor  all  ions  since  1  ~  v.  cannot  be  constant  as  a  function  of  r  when  transverse  velocity 
components  exist  Consideration  is  therefore  limited  here  to  cases  where  I  '  ^  /  ,c  cascs 
for  which  the  transverse  energy  content  of  the  beam  is  small 
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The  constants  A,  B ,  C,  and  <f>  are  chosen  to  reflect  the  desired  velocity  distribution  at  z 


0. 


AE  +  Bp  +  CU(z  =  0)  =  (vr  +  rVjif )2  +  ( v„  -  cur)2 
For  an  initially  diverging  beam,  z,  is  negative.  Equation  (6)  requires  that 
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where  <r  =  Vjilz The  small-transverse-energy  constraint  requires  that  consideration  be 
limited  to  beams  injected  at  small  angles  to  the  guide  field,  i.e. 

|z,|  »  R  ,  ilR  «  V„ 
for  the  radial  and  azimuthal  velocity  components. 

The  radial-density  distribution  at  z  =  0  can  now  be  determined  by  integrating  / over  velo¬ 
city  space.  This  exercise  results  in  the  expression 
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and  S  =  RU/y,„.  Density  profiles  for  various  values  of  A  are  shown  in  Fig.  (2).  For  the 
"cool"  beams  of  interest  here  (  F„  »  F,„),  the  beam  profile  is  nearly  square. 


The  average  beam  area  as  a  function  of  axial  position  can  be  determined  by  integrating  / 
evaluated  along  the  axis  of  symmetry.  This  calculation  yields 
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where 
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Equations  (9)-(U)  show  that  for  an  initially  diverging  beam,  n(0,;)  first  decreases  from 
n„  and  then  increase  as  the  beam  approaches  a  focus.  The  focus  location  L ,  the  value  of  ;  for 
maximum  density  occurs  very  close  to  the  point  where  a  is  minimum.  From  Eq.  (11), 
minimum  a  occurs  at  5  —  5,.  where 


I'-aL/V'  +  i-lZ- 


with  6  determined  from  the  value  of  it  and  Eq.  (7).  The  radial  beam  density  distribution  at 
5  ”  5,  can  again  be  determined  by  integration  of /over  velocity  space  with  the  result 


n(r.L) 


A  exp  —  A- 


n(o,L) 


is  obtained  from  Eq.  (9)  with  5  -  5,  The  quantity  A  is  a  measure  of  the  degree  of  current- 
density  multiplication  occuring  at  the  best  focus.  The  variation  of  A  with  the  thermal  parame¬ 
ter  8  is  shown  in  Fig  3  lor  various  values  of  ir  and  the  two  limiting  rotational  frequencies. 
Best  focussing  is  associated  with  to  -  -  Vi  ft  since  this  case  corresponds  to  zero  canonical  angu¬ 
lar  momentum.  For  this  choice,  the  minimum  value  of  0(5)  is  zero  so  that  A  may  be  obtained 


from  Eq.  (9)  in  the  limit  of  a  going  to  zero. 


A\to  — 


4M 


54 


MOSHER  el  at. 


As  suggested  earlier,  this  is  the  same  degree  of  focussing  one  would  obtain  ballistically  using 
the  focal  length  z,  rather  than  L  in  Eq.  (1). 

As  an  example  of  how  reflexive  focussing  can  be  employed  in  a  pellet-fusion  reactor,  con¬ 
sider  a  10  MA  beam  of  10  MeV  protons  focussed  onto  a  pellet  of  1  cm2  cross-sectional  area  5  x 
102  cm  away  from  the  diode.  The  required  degree  of  focussing  is  A  —  nR2  where  R  is  the 
ion-emitting  radius  of  the  anode.  Assuming  that  the  m^jor  contribution  to  V,h  arises  from 
fluctuating  fields  in  the  diode  leading  to  an  about  ±1°  time-dependent  deflection  of  the  beam- 
emission  angle,  then  VjV,h  =  60  so  that  Eq.  (15)  requires  z  =  30  cm.  Simultaneous  solution 
of  Eqs.  (7)  and  (12)  leads  to  tlLIV0~  5  so  that  the  required  solenoidal  field  is  about  5  kG. 
Choosing  R  =  5  cm  in  order  to  satisfy  the  small-transverse- velocity  requirement  |zA|  »  R 
leads  to  a  diode  aspect  ratio  of  15  for  a  3  mm  diode  spacing.  This  ratio  is  well  matched  to  the 
generator  required  to  genera'  the  1  ft  ion  beam. 
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V.  MHD  RESPONSE  OF  TRANSPORT  CHANNELS 

A.  INTRODUCTION 

With  recent  advances  in  the  technology  of  intense  light-ion  beam  generators,1  2  the 
suggestion  that  these  ion  beam  could  be  used  in  a  pellet  fusion  scheme  has  drawn  considerable 
interest3  as  have  similar  concepts  using  heavy  ion  beams4  (e  g.  20  GeV  U+  ions).  In  a  func¬ 
tioning  reactor,  the  ion  beam  would  have  to  be  propagated  a  distance  L  (about  10  m)  from  the 
acceleration  region  while  employing  some  focussing  scheme  to  deliver  the  beam  to  the  target. 

For  heavy  ion  beams,  field-free  ballistic  focussing  to  the  target  has  been  proposed.4  5  For 
light  ion  beams  (e.g.  10  MeV  protons),  one  concept6  that  is  under  investigation  involves 
focussing  the  beam  just  outside  the  diode  region  by  known  techniques'  and  then  injecting  the 
focussed  beam  into  a  z-pinch  plasma  channel.  The  beam  remains  focussed  during  propagation 
by  the  magnetic  field  in  the  plasma  channel  and  is  delivered  on  target  with  no  additional  focuss¬ 
ing.  For  both  of  these  techniques,  it  must  be  demonstrated  that  the  self-consistant  electromag¬ 
netic  fields  associated  with  the  plasma  response  to  beam  passage  do  not  disrupt  transport  and 
focussing.  We  have  started  to  investigate  this  problem  by  concentrating  on  the  physics  occur¬ 
ring  in  the  plasma  channel  as  a  result  of  the  ion  beam  passage.  This  area  of  investigation  has 
been  largely  ignored  in  studies  so  far;  it  has  been  generally  assumed  that  the  channel  does  not 
respond  to  the  beam  on  the  time  scale  involved.  With  beam  pulse  lengths  of  50  ns  for  light  ion 
beams  and  of  around  10  ns  or  more  for  heavy  ion  beams,  this  assumption  is  not  generally 
correct  for  transport  in  channels  of  cm  radius  or  less.  Magnetic  nonneutralization  effects  are 
important  for  large  radius  beams  as  well,  since  small  (i.e.  —  1  kA)  net  currents  are  sufficient  to 
destroy  focussibility  to  mm  spot  sizes  for  heavy  ion  beams  which  are  ballistically  focussed. 
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Recent  work  has  concentrated  on  the  plasma  response  to  light-ion  beams  with  a  view  to 
answering  the  questions: 

(a)  Does  a  plasma-parameter  window  exist  which  allows  efficient  transport  over  10  meter 
distances? 

(b)  Do  the  self-consistent  electromagnetic  fields  associated  with  beam  heating,  net  current 
effects  and  plasma  hydrodynamics  disrupt  beam  transport  or  its  ability  to  bunch? 

Some  constraints  which  have  to  be  satisfied  by  the  beam  plasma  system  include: 

(a)  The  density  in  the  hole  bored  by  the  applied  z-discharge  and  beam  currents  must  be 
sufficiently  high  to  permit  complete  charge-neutralization. 

(b)  The  beam-heated  plasma  must  be  sufficiently  high  in  temperature  to  provide  good 
magnetic  neutralization. 

(c)  The  temperature  must  be  sufficiently  high  to  reduce  return-current  electric  fields  to  a 
level  which  neither  slows  the  beam  appreciably  nor  disrupts  bunching. 

B.  SUMMARY  OF  CODE  RESULTS 

A  1-D  radial,  time-dependent,  Eulerian  code  has  been  developed  recently  to  answer  these 
questions.  The  code  solves  the  plasma  moment  equations.  Maxwell  equations,  and  transport 
equations  using  coefficients  for  electrical  and  thermal  conduction  taken  from  Braginskii.7  The 
beam  is  specified  by  its  current  density  in  space  and  time.  Direct  beam  energy  deposition  and 
radiation  cooling  as  well  as  return-current  heating  have  been  included  in  the  code.  Initial  con¬ 
ditions  consist  of  a  plasma  z-discharge  in  a  high-density  gas.8  The  density  of  the  discharge 
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channel  is  typically  .01  -.1  atmosphere  of  low-z  gas  while  the  background  density  may  be  much 
higher. 

Preliminary  results  appear  in  Fig.  1.  They  correspond  to  a  1  MA,  10  MeV  proton  beam, 
0.6  cm  in  radius  propagating  in  a  preformed  hydrogen  channel  of  the  same  dimension  with  an 
initial  density  on  the  axis  of  2  x  1018.  Results  shown  are  radial  prohles  for  density,  azimuthal- 
magnetic  and  axial-electric  fields,  and  temperature  at  t  —  40  ns  into  the  beam  pulse.  These  are 

the  conditions  which  develop  near  the  tail  of  the  beam.  No  impurities  have  been  included  in 
this  calculation  so  that  only  hydrogen  bremsstrahlung  losses  are  included.  Time  histories 
through  any  point  fixed  in  the  chamber  of  the  maximum  B-field,  maximum  temperature,  cen¬ 
tral  electric-field  strength  and  minimum  channel  density  appear  in  Fig.  2. 

Conclusions  which  can  be  drawn  at  this  point  are: 

(a)  Hydrogen  backgrounds  of  1-3  x  10l8cm~ 3  allow  propagation  satisfying  all  constraints 
to  transport  and  bunching. 

(b)  Hydrodynamic  expansion  at  lower  densities  may  result  in  insufficient  neutralizing  den¬ 
sities. 


(c)  Low  plasma  temperatures  at  higher  densities  cause  excessive  return-current  electric 
field  strengths  which  slow  the  beam  and  perturb  bunching.  Line  radiation  from  impurities  in 
the  background  may  also  cause  low  temperatures  to  develop  if  present  in  sufficient  densities. 

Phenomena  which  have  not  been  included  in  the  code  yet  and  which  are  important  to  a 
more  complete  understanding  of  the  channel  physics  are: 

•  Non-equilibrium  ionization  of  the  plasma  at  the  head  of  the  beam. 
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•  Values  for  electrical  conductivity  at  low  temperature  which  determines  the  penetration  of 
the  B-field  into  the  plasma  and  the  degree  of  current  non-neutralization  achieved  late  on 
in  the  pulse. 

These  phenomena  plus  comparison  with  experiments  at  lower  current  intensities  (with 
correspondingly  lower  temperatures  and  ionization  levels),  extension  of  this  work  to  heavy  ion 
beams,  and  a  more  consistent  treatment  of  the  ion-current  distribution  in  position  and  time  due 
to  self  fields  and  bunching  will  constitute  future  developments  of  this  work. 
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